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The symmetries of the DNA double helix require a new term in its linear response 
to stress: the coupling between twist and stretch. Recent experiments with torsionally- 
constrained single molecules give the first direct measurement of this important material 
parameter. We extract its value from a recent experiment of Strick et al. [Science 271 
(1996) 1835] and find rough agreement with an independent experimental estimate recently 
given by Marko. We also present a very simple microscopic theory predicting a value 
comparable to the one observed. 
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1. Introduction 



The idea of studying the response of DNA to mechanical stress is as old as the discov- 
ery of the double helix structure itself QTj] . While many elements of DNA function require 
detailed understanding of specific chemical bonds (for example the binding of small lig- 
ands), still others are quite nonspecific and reflect overall mechanical properties. Moreover, 
since the helix repeat distance of £ ~ 3.4 nm involves dozens of atoms, it is reasonable to 
hope that this length-scale regime would be long enough so that the cooperative response 
of many atoms would justify the use of a continuum, classical theory, yet short enough that 
the spatial structure of DNA matters. In this Letter we will argue that this expectation is 
indeed fulfilled. 

Since moreover various important biological processes involve length scales comparable 
to £ (notably the winding of DNA onto histones), the details of this elasticity theory 
should prove important. Yet until recently little was known about the relevant elastic 
constants. Extensive experimental work yielded fair agreement on the values of the bend 
and twist persistence lengths, though the former was plagued with uncertainties due to 
the polyelectrolyte character of DNA 0. A simple model of DNA as a circular elastic 
rod gives a reasonable account of many features of its long-scale behavior, for example 
supercoiling 0[Q. Some authors have sought to justify this model by invoking a shell of 
structured water around the DNA ||. 

Recently, techniques of micromanipulation via optical tweezers and magnetic beads 
have yielded improved values for the bend stiffness from the phenomenon of thermally- 
induced entropic elasticity || , as well as the direct measurement of a third elastic con- 
stant, the stretch modulus, by exploring the force range 10-50pN J7|[|5J. Significantly, the 
relation between bending stiffness, stretch modulus, and the diameter of DNA turned out 
to be roughly as predicted from the classical theory of beam elasticity @ |§ |§ , supporting 
the expectations mentioned above. 

Still missing, however, has been any direct measurement of the elastic constants reflect- 
ing the chiral (i.e. helical) character of DNA. One such constant, a twist-bend coupling, 



was investigated by Marko and Siggia ||10|1 , but no direct experimental measurement has 
yet been devised. In this Letter we introduce a new chiral coupling, the twist-stretch en- 
ergy. Electrostatic effects do not complicate the analysis of this coupling. We will explain 
why our term is needed, extract its value from the experiment of Strick et al. ||11|| , and 



compare it to a the prediction of a simple microscopic model to see that its magnitude 

1 



is in line with the expectations of classical elasticity theory. J. Marko has independently 
introduced the same coupling and estimated its value from different experiments [|T2f : 
values are in rough agreement. 
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2. Experiment 

DNA differs from simpler polymers in that it can resist twisting, but it is not easy to 
measure this effect directly due to the difficulty of applying external torques to a single 
molecule. Early investigations of DNA twist were either limited to passive, fluorescence- 
depolarization measurements @ , or else to studying global shape changes in circular DNA 
of varying linkage ||. The first single- molecule stretching experiments constrained only 
the locations of the two ends of the DNA strand. The unique feature of the experiment of 
Strick et al. was the added ability to constrain the orientation of each end of the molecule. 

We will study Fig. 3a of ref. ||11|| . In this experiment, a constant force of 8pN was 
applied to the molecule and the end-to-end length z tot monitored as the terminal end was 
rotated through ALk turns from its relaxed state (which has Lk turns). In this way the 
helix could be over- or undertwisted by as much as ±10%. Over this range of imposed 
linkage z tot was found to be a linear function of a: 

e = const. — 0.15a where a = ALk/Lk and e = (z tot /z tot ) — 1 . (2-1) 

Thus a is the fractional excess link and e is the extension relative to the relaxed state. 
Eqn. (|2.1] ) is the experimentally observed twist-stretch coupling. 

The existence of a linear term in ( |2.1|) is direct evidence of the chiral character of 
the molecule, and its sign is as expected on geometrical grounds: untwisting the molecule 
tends to lengthen it. Still geometry alone cannot explain this result. Consider the outer 
sugar-phosphate backbones of the DNA. Suppose that the twist-stretch phenomenon were 
due to the straightening of these helical backbones while they maintained constant length, 
0.6 nm per phosphate, and constant distance 0.9 nm from the center of the molecule. 
Then since each basepair step is h = 0.34 nm high, the circumferential length per step 



is £ c = V-6 2 — -34 2 nm. The corresponding twist angle per step is given by [|13] 9 = 
(£ c /2)/.9nm = 32°, roughly as observed. Supposing now an extension by Ah/h = e, we 
find an untwisting by a = 69/9 = const. — e/2.0, quite different from what is observed, 
eqn. ( |2.1| ). We must seek an explanation of the experimental result not in terms of a 
geometrical ball-and-stick model but in the context of an elastic response theory. 
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3. Simple Model 

We will begin by neglecting bend fluctuations (see below). A straight rod under 
tension and torque will stretch and twist. We can describe it by the reduced elastic free 
energy 

fM e) = = ^ [Ca 2 + Be 2 + 2Dea] - re . (3.1) 

The twist persistence length is C ~ 75 nm ||, while the helix parameter uj = 2tt/£ = 
1.85/nm. We will takeB B m 1100 pN/uJ 2 k B T m 78 nm |J. In the experiment under study 
the reduced force is r = 8pN//c s T « 1.95/nm. For a circular beam made of isotropic 
material the cross-term D is absent || , since twisting is odd under spatial inversion while 
stretching is even. For a helical beam, however, we must expect to find this term. 

We now minimize f x with respect to e, a at fixed tension with an imposed constraint 
on the overtwist a. Minimizing at fixed a and r gives A = uj 2 (Ds + Ca) and hence 

e = e a=0 - (D/B)a . (3.2) 



Comparing to (|2.1| ), we obtain the desired result: D = 12 nm. To compare this to Marko's 
analysis, we note that his dimensionless g equals our Du , so that we get g = 22. The 
rough agreement with Marko's result g = 35 [0 indicates that the data show a real 
material parameter of DNA and not some artifact. We do not expect exact agreement, 
since Marko studied the nonlinear overstretching transition of |14|] ; our value came from 
the linear regime of small strains. 



4. Bend Fluctuations 

To arrive at ([O]) we listed the variables which were constrained, coupled to external 
forces, and/or observed in the experiment, namely e and a, then wrote the most general 
quadratic function allowed by symmetry. Thus (|3.1| ) is a phenomenological model; its 
coefficients C, B, D reflect both intrinsic elasticity and the effects of thermal fluctuations. 

1 B reflects the intrinsic stretchiness of DNA, since electrostatic self-repulsion simply 
shifts the equilibrium length without affecting the spring constant. Indeed experiments 
show little or no dependence of B on salt, unlike the situation with the effective bend 
persistence length |§. We also expect C to reflect intrinsic elasticity, since twisting does 
not affect the average charge distribution. 
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Indeed it is well known that thermal bend fluctuations reduce the effective stretch modulus 
at modest tension via the "entropic elasticity" effect ||||. Somewhat inconsistently we 
arrived at our value of D by using the intrinsic stretch modulus in (J3.2|). In this section we 



will justify the procedure by introducing a more elaborate model with bend fluctuations, 
rederiving the analog of ( |3.2| ), and again comparing to ( |2.1| ). 

We begin by defining local variables (Fig. la) (see |T^[]Tt|). DNA is a stack of base 
pairs. We neglect sequence effects and so regard all base pairs as copies of one standard 
slice. The standard slice contains a reference point with the property that the locus of these 
is the helix axis, a straight line of length L in the relaxed state. Through this reference 
point we next draw a fixed vector; a convenient choice is the "dyad" pointing into the 
minor groove and perpendicular to the helix axis. 



(b) 

Fig. 1: Schematic diagrams defining variables used in the text. The 
offset from the helix axis has been exaggerated for clarity, a) Notation 
used in the study of bend fluctuations. We describe the DNA by the 
helix axis (dotted curve) and the axis E x , which is a fixed vector in each 
base pair, b) Notation used in the microscopic model. The helix axis 
(dotted line) is now supposed straight. We describe the DNA by the 
dashed curve and the axis E 1 as before. 
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To describe stressed states, we simply specify the locus of reference points as a pa- 
rameterized curve in space (dotted line in Fig. la) and the dyad as a field of vectors E 1 
normal to this curve. We let E 3 be the unit tangent to the axis and complete to an or- 
thonormal triad by defining E 2 = E 3 x E 1 . Next we introduce a parameter s to label each 
slice; s corresponds to arc length along the original, unstressed helix axis and so always 
runs from to L. The actual arc length along the distorted axis will not be ds but rather 
(1 + a(s))ds; a is thus the axial strain. 

Thus our local variables are E^s) and a(s). Our program consists of four steps: 
i) Find the strains in terms of the local variables, ii) Write the general linear elasticity 
theory of these strains with a force coupling to the extension e and a torque coupling to 
the twist a. Hi) Compute and minimize the free energy to find the end-to-end length (z tot ) 
in terms of the constrained a and the applied force r. (It will prove convenient in this step 
to convert to new variables t^s), t 2 (s), <p(s), and a(s).) iv) We will then be able to relate 
the experimental result to intrinsic elastic constants. We will see that our naive calculation 
of the previous section is justified. Details of this calculation will appear elsewhere. 

Step i: In the relaxed state each slice bears a constant relation to its predecessor. 
Thus while E 10 , E 20 , and E 30 all vary in space, the derivatives with respect to s are of 
the form 

dEjds = -e^ k n j0 E k0 , (4.1) 

where O - are the constants (0,0, u> ). For the deformed state, (|4.1|) defines the functions 
f^(s). Our strain variables are then 1 (s), 2 (s), 3 (s) — u , and a(s). 

Step ii: Our strain variables have the desirable property that in terms of them the 
elastic constants are independent of s. Moreover the end-for-end symmetry of DNA implies 
that the elastic matrix is unchanged upon changing the sign of 1 [|IU|. Thus we generalize 



the model of JlO] to 



I r L r 

f 2 = ds A'Slj 2 + AQ 2 2 + C(0 3 - uj ) 2 + Buj 2 a 2 + 2L>w (fi 3 - uj Q )a 

2L Jo i - ' (42) 

+ 2G(0 3 - u) )tt 2 + 2Kuo o n 2 a . 



Here G is the twist-bend coupling of | 10[j , while K is an allowed coupling between stretch 
and bend. 

We will apply a perturbative treatment to (|4.2| ) . Such an approximation is valid since 
in the experiment we are analyzing the applied force is large enough to keep the end-to- 
end distance over 90% the full relaxed contour length, but not large enough to create large 
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intrinsic stretch a. In addition, the applied overtwist a is at most ±10% \'\ I]; indeed the 
slope reported in ( |2.1| ) can be found from an even smaller range of a than this. 

Step Hi: It proves useful to refer the frame {E 1 , E 2 , E 3 } to a reference frame 
{e l5 e 2 , z} rotating at u : 



e x = xcos(u s) + y sin(w s) e 2 = —xsm(u s) +ycos(u> s) . 



(4.3) 



We will then write the deformed frame in terms of three small quantities: two deviations 
of the tangent vector t 1 2 and an angle (p. To first order in these we have 



E 1 = e x + fe 2 — t\Z , E 2 = —fe 1 + e 2 — t 2 z , E 3 = t 1 e l + t 2 e 2 + z 



(4.4) 



Thus E 1 is the unit vector perpendicular to the tangent whose projection to the xy plane 
makes angle <p with the reference frame. The advantage of these coordinates is now that 
the excess link is simply the difference <p(L) — <p(0). The relation between Link, Twist, 
and Writhe will also emerge automatically instead of having to be enforced by hand. 

Exercising the definitions we find that in terms of £ l5 t 2 , ip the strains are 1 = 
—i 2 — u>) Q t x , fl 2 =t 1 — uJ t 2 , fi 3 = iv + if. Substituting into ( |4.2[ ), adding external tension 
r and torque A, and rearranging gives the free energy 



f2 = 1 kJ ds A '( i 2+ UJ o t i) 2 + A{i 1 -u t 2 + jf+^ay 



+ C 



91 
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(4.5) 



T - J d S (l + a)(l-it 1 2 -it 2 2 ) J dsy. 



Step iv: Let us first focus only on the linear couplings to the applied forces r, A. 
Then eqn. ( f4.5|) shows that the Lagrange multiplier A couples only to the constant part of 
if, which is just the overtwist u a. Also the constant part of a is the extension e. Thus 
neglecting the quadratic couplings to r gives that only the constant modes t l5 t 2 , cr, and 
e respond to r and cr, and they minimize the function 



iOr 



A'tf + A[t 2 
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(4.6) 



We see that t 2 , which is neither controlled nor observed, adjusts slightly but that the 
measured twist-stretch coupling D reflects the combination D — ^jf- of intrinsic elastic 
parameters. Parenthetically we note that C = C — is nearly equal to C because we 
expect G to be small, and similarly for B. The corrections are small because they reflect 
the small deviation of DNA from a straight circular rod. 

We are almost done. Normally the linear force terms suffice since the applied force 
t/uj q 2 is small compared to the persistence length A. The only exception to this rule 
comes from those modes which decouple completely from ( |4.5| ) when r = 0: these modes 
will have large fluctuations, and the nonlinear couplings to r are needed to cut them off. 
These dangerous modes are the Fourier modes of t 1 , t 2 of wavenumber close to u , as we 
see from ( |4.5|) . Their thermal fluctuations indeed make a large contribution to (z tot ) at 
small force. Fortunately, though, this contribution is completely independent of a, since a 
enters linearly in the problem and hence makes no contribution to the fluctuation problem. 
Thus we can simply absorb this contribution to (z tot ) into the constant term of (|2.1|) . 

Thus we have found the interpretation of the experimentally-determined twist-stretch 



coupling found in section 2: in terms of the intrinsic elasticity of DNA the slope in ( |2 . 1| ) 
fixes the combination (DA — GK) /AB to be 22 nm. The low-force stretching experiments 
give bend stiflhessi A = 40 nm ||. Other experiments fix 5, C to the approximate values 
quoted earlier. The remaining two combinations of the six couplings in ( f4.5|) do not appear 
to be relevant for existing experiments. 



5. Microscopic Model 

The elastic theory in the previous section was very general, but it gave no indication 
of the expected magnitudes of the various couplings. To gain further confidence in our 
result, we will now see how to estimate the expected twist-stretch coupling based on the 
measured values of the other elastic constants and geometrical information about DNA. 
We will use a simple, intuitive microscopic picture of DNA as a helical rod to show how 
twist-stretch coupling can arise and get its general scaling with the geometric parameters. 
While the model is unrealistic it captures the underlying symmetries and shows that the 
value of D calculated above is reasonable. 

2 This value for A is slightly smaller than the traditional one. The authors of ref. || 
eliminated the electrostatic contribution to A by extrapolating to high concentration of 
high- valence added salt. 
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Our picture will be a beam of isotropic elastic material of circular cross-section, ini- 
tially bent into a helix of pitch £ , with the beam center slightly displaced from the helix 
axis by d <C £ Fig. 2b defines notation (for realistic depictions of DNA structure see 
for example [|13|])- In this section we will consider only uniform deformations of the helix; 
in particular the helix axis will always be straight. It proves convenient to define slightly 
different variables from the previous section: instead of following the helix axis, now our 
curve follows the centerline of the beam. We again call the tangent to this curve E 30 (s), 
where the arc length s runs from to L, but now L is slightly longer than the end-to-end 
length of the relaxed beam. Next we draw a second curve, the locus of points farthest from 
the helix axis. Let E 10 (s) be the field of vectors perpendicular to the tangent E 30 (s) and 
pointing from the first to the second curve. Finally complete E 30 , E 10 to an orthonormal 
triad by defining E 20 = E 30 x E 1Q . 

The distorted beam will then have its modified frame {E^s), E 2 (s), E 3 (s)}, where 
now s is the arc length rescaled by (1 + a) -1 to again run from to L and a is the axial 
strain as before. We also define strain variables as before; for the uniform deformations 
considered these are constants independent of s. For a nearly-straight, circular rod the 
elastic energy is then M 

f 3 = \ [A(Q 2 - O 20 ) 2 + C(n 3 - O 30 ) 2 + Bco 2 a 2 ] . (5.1) 

We have rotated our reference frames about the tangent vector to eliminate Q x . The 
constants A, C, and B can in turn be expressed in terms of the effective Young modulus, 
shear modulus, and diameter of the rod, but we instead use the measured values quoted 
earlier. 

To use ( |5.1|) we need to find fi 2 an d ^3 in terms of the helix parameters: helix 
axis offset <i, end-to-end length z tot , and total rotation of the cross-section. The latter 
quantity plays the role of linking number for open DNA, and so we will call it Lk. To 
get the required relations it is helpful to use the physical image of a gyroscope rotating 
at "angular frequency" about an axis parallel to O while moving at constant "speed" 
along an axis E 3 fixed in the body. We take "time" to run from to L, the original relaxed 
arc length; to allow for intrinsic stretching we then take the "speed" to be 1 + a. One then 
finds that 

d = n 2 (l + a)/\n\ 2 , z tot = Zn 3 (l + a)/\n\, Lk = L|fi|/27r. (5.2) 
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To fix O we impose the values a = 0, u> = 27rLk / z tot = 1.85/nm, and a helix offset 
d . We will choose d to get the observed value of D and see that it is reasonable. Working 
to second order in d ( |5T2"1 ) gives O 30 = u Q (l — c^Wq 2 ), O 20 = <i c<; 
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We must now minimize with the constraint of fixed z tot and Lk. Let z tot = 

(1 + e)z tot , so that e again measures changes in end-to-end distance, and Lk = (1 + a)Lk 



as usual. Again using (|5.2|) , one finds 



^2-^20=^0' fi 3 - fi 30 = u a ~ ^oX/ 3 > « = £ + u O d oP ~ ( w O rf o) 2cr > ( 5 - 3 ) 



where (3 is free. Substituting into ( |5.1|) and minimizing over (3 reveals a as coupling 
corresponding to D = (u d ) 2 (C — A)B/A. This formula fits our measured value of D if 
we choose d Q = 0.2 nm. 

The value of d Q is not known a priori, since of course DNA is not really an elastic con- 
tinuum with circular cross-section. Nevertheless, inspection of the known molecular struc- 
ture indeed suggests an elastic center offset from the helix axis by a couple of Angstroms 
||13|| . In any case we have shown that the measured value of D is of the order of magnitude 
expected from simple elasticity theory.i 



6. Conclusion 

We have pointed out a strong twist-stretch coupling in torsionally-constrained DNA 
stretching experiments, evaluated it, argued that it reflects intrinsic elasticity of the DNA 
duplex, and shown that the value we obtained is consistent with elementary considerations 
from classical elasticity theory. A greater challenge remains to predict this coupling from 
the wealth of available crystallographic information on the conformation of short oligomers. 
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3 Actually a helical beam can have a twist-stretch coupling even if its axis is on center, 
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